High T c superconductivity due to coexisting wide and narrow bands: 
A fluctuation exchange study on the Hubbard ladder as a test case 
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We propose that when the Fermi level lies within a wide band and also lies close to but not within 
a coexisting narrow band, high T c superconductivity may take place due to the large number of 
interband pair scattering channels and the small renormalization of the quasiparticles. We show 
using fluctuation exchange method that this mechanism works for the Hubbard model on a ladder 
lattice with diagonal hoppings. From this viewpoint, we give a possible explanation for the low T c 
for the actual hole doped ladder compound, and further predict a higher T c for the case of electron 
doping. 
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The discovery of high T c superconductivity in the 
cuprates 1 , followed by discoveries of various unconven- 
tional superconductors, has brought up renewed fascina- 
tion for the search of high T c superconductors, for which 
a theoretical guiding principle is highly desired. One way 
to attack this problem is to theoretically search for lat- 
tice structures that provide good conditions for Cooper 
pairing. Ladder-like structures may be considered as a 
candidate for this^, but up to now, T c in the actual hole 
doped ladder compound Sri4_ x Ca x Cii2404i (14-24-41) 
remains to be around ~ 12 K, 3 which is low compared to 
the layered cuprates. In this context, we have previously 
proposed a high T c mechanism due to disconnected Fermi 
surfaces, which can in fact be realized with a ladder-like 
lattice structure with a larger hopping integral in the 
rung direction than in the leg direction^ but up to now 
there exist no actual materials (or methods) to realize 
the situation we have proposed. In the present study, we 
propose a different mechanism for high T c superconduc- 
tivity in systems with coexisting narrow and wide bands, 
which may be realized in actual ladder compounds with 
electron doping. 

Our idea is as follows. Let us consider a system where 
the Fermi level Ep lies within a band with a moderate 
width (we call this "the wide band A" ) and also lies close 
to, but not within, a narrow band ("band B"). If the am- 
plitude of the pair scattering processes from band A to 
band B is strong enough, the pairing instability may be- 
come very large because (i) the sign change of the gap, 
necessary for repulsive pairing interactions, occurs be- 
tween the bands, so that there are no gap nodes on the 
Fermi surface, and (ii) there is a huge number of inter- 
band pair scattering channels due to the narrow character 
of band B (see the bottom of the right panel in Fig^l . 
Another point to be stressed here is that the Fermi level 
is not within the narrow band, so that the renormaliza- 
tion of the quasiparticles at the Fermi level is not so large 
as to strongly suppress superconductivity. Systems with 
coexisting wide and narrow bands may be reminiscent of 
models consisting of wide s or p bands and narrow d or 
/ bands, but in that case, the amplitude of the inter- 
band (interorbital) pair scattering processes may not be 




FIG. 1: Left panel: The hopping integrals in the ladder 
lattice is shown. Right panel : the energy dispersion for the 
ladder lattice with t' = (top) and with t' — —ti (bottom). 
The curved arrows show the pair scattering processes that 
give rise to superconductivity. 



so large because of the different character of the orbitals. 
Here we consider systems with one orbital per site, where 
the multiplicity of the bands originates from the lattice 
structure rather than the multiplicity of the orbitals at 
each site, so that the pairing interaction can be strong, 
originating from the large on-site repulsion. 

The above condition for the energy bands can most 
simply be satisfied in a tightbinding model on a lad- 
der lattice with diagonal hoppings shown in Fig^ The 
Hamiltonian of this model is given in momentum space 
as 



kin 



= £(4 



ka ' a ka 



ka 



—2ti cos k 
-2t' cos k — t r 



-2t' cos k — t r 
—2ti cos k 



Cka 
dker 



where Cka and d^ a annihilate an electron with spin a 
at wave number k on the left and right legs, respectively, 
and ti,t r , and t' are the hopping integrals in the leg, rung, 
and diagonal directions, respectively. The dispersion of 
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FIG. 2: The trellis lattice adopted in the present study. 



the two bands are given as s±(k) = —2(ti ± t') cos k =F t r . 
When t' = (FigQ top of the right panel) , the two bands 
have identical dispersions with a level offset of 2t r , while 
one of the bands (band B) is narrower than the other 
(band A) in the presence of t' , and becomes perfectly 
flat for t' = ±t ; (Fig^ bottom of the right panel) . 

Here, we consider the on-site interaction (U) term in 
addition to the above kinetic energy terms, and also take 
into account the trellis-like lattice structure of the actual 
cuprate ladder compounds, where the ladders are weakly 
coupled by diagonal hoppings We estimate the su- 
perconducting transition temperature of this Hubbard 
model using the combination of the fluctuation exchange 
method (FLEX) and the Eliashberg equation, which has 
been successfully applied to the problem of layered high 
T c cupratesiiAii The values of U, t r , t% will be fixed at 
U = 6ti, t r = ti, and ij = 0.25t; throughout the study. 
We define the the band filling n as n — [number of elec- 
trons] /[number of sites], so when the bands are both fully 
filled, the band filling is n = 2. 

In the two band version of FLEX^ the Green's func- 
tion G, the susceptibility %, the self-energy E, and the 
superconducting gap function all become 2x2 matri- 
ces, e.g., G; m (fc, ie n ), where I, m specify the two sites in a 
unit cell. The orbital-indexed matrices for Green's func- 
tion and the gap functions can be converted into band- 
indexed ones with a unitary transformation. As for the 
spin susceptibility, we diagonalize the spin susceptibility 
matrix and concentrate on the larger eigenvalue, denoted 
as x- 

The actual calculation proceeds as (i) Dyson's equa- 
tion is solved to obtain the renormalized Green's func- 
tion G(k), where k = (k,ie n ) denotes the 2D wave- 
vectors and the Matsubara frequencies, (ii) the effec- 
tive electron-electron interaction V^\q) is calculated by 
collecting RPA-type bubbles and ladder diagrams con- 
sisting of the renormalized Green's function, namely, 
by summing up powers of the irreducible susceptibil- 
ity xirr(g) = -j? E fc G(fc + q)G{k) (^number of fc- 
point meshes), (iii) the self energy is obtained as = 
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FIG. 3: T c plotted (a) as a function of — t' for n — 1.25, and 
(b) as a function of n for — t 1 = 0.95t;. (c)The eigenvalue 
of the Eliashberg equation A plotted as a function of —t' for 
n = 0.9 and T = 0.05U. 



J2 q G(k — q)V^(q), which is substituted into Dyson's 
equation in (i), and these procedures are repeated until 
convergence is attained. 

We determine T c as the temperature at which the 
eigenvalue A of the linearized Eliashberg equation, 

\cj> lm (k) = -^ J2 ^^-k')Giv{k')G mml {-k')cj> llm ,{k'), 

h'l'm' 

reaches unity. Here the pairing interaction for sin- 
glet pairing is given by V^(q) = U + '^U 2 Xm(q)/(l — 

UXirr(q)) ~ \U 2 XuM)/{^ + UXuM))- 

Throughout the study, we take up to 64 x 64 fc-point 
meshes and the Matsubara frequencies e„ from — (2N C — 
1)ttT to (2N C - 1)ttT with iV c up to 8192 in order to 
ensure convergence at low temperatures. 

We now move on to the results. We first show the t' 
dependence of T c for n — 1.25. We consider the case 
of t' < since this is the realistic choice of sign for the 
cupratesiiS In Fig|3Ja), we plot T c as a function of — t' 
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X(kx,k y ,0), (c) (j> a (k x ,k y ,ink B T) for n= 1.25, t 



(b) 
0.95t h 



and T = 0.08t; viewed from the direction of the k x axis. The 
thickness of the curves represents the dispersion in the k x di- 
rection. In (a) and (c), quantities for a = u (upper band) and 
I (lower band) are shown by dark green and red curves, respec- 
tively, while N.B (W.B.) denotes the portions corresponding 
to those for the narrow band B (wide band A). The relation 
between the narrow/wide bands and the lower/upper bands 
is shown in the inset of figure (a). 



for n — 1.25. It can be seen that T c takes its maximum 
around —if = ti, where band B is flat. There, T c almost 
reaches 0.08*;, which is extremely high if ti is assumed to 
be of the order of few hundred meV as in the cuprates. 

To trace back the origin of this high T c , we look into the 
Green's functions, spin susceptibility and the gap func- 
tions for —t' = 0.95*; and n = 1.25. In Fig0] we plot 
IGujikx^kyJirT)^, x( k x,k y ,0), and </) Ut i(k x ,k v ,iirT) 
viewed from the direction of the k x axis. Here, we display 
the results for G and <fi in the band representation, where 
u and I denote the upper and lower portions of the bands 
as shown in the inset of Fig^a). Since bands A and 
B intersect with each other, the wide band A (the nar- 
row band B) consists of the upper (lower) band around 
k y = ±7r and the lower (upper) band around k y = 0. In 
the Green's functions, the two peak-like structures seen 
around k y — ±7r/4 is due to the Fermi level crossing, 
while the narrow structure (noted as N.B.) owes to the 
flatness of band B. Since the volume of the Fermi surface 
(the length 2fcp) is ~ 7t/2, the wide band A is nearly 



IGI 2 o.( 



(a) 




















- 
- 










-^i i i^- 


-71 


-It/2 




ky 


Till 


Jt 




(c) 


1 -J 







-ji/2 



71/2 



FIG. 5: Plots similar to FigEJfor n = 0.9 and t' = -0.95t;. 



0.06 
0.05 
0.04 
T c /t, 0.03 
0.02 
0.01 
0.00 



/ 




1.1 1.2 

n 



FIG. 6: T c plotted as a function of n for t' 



-0.4ti 



quarter filled, and thus the narrow band is fully filled, 
i.e., the Fermi level lies above the narrow band. The 
spin susceptibility has a broad structure again due to the 
flatness of band B, which enhances the number of pair 
scattering channels due to spin fluctuations. The gap 
function changes sign between band A and band B, but 
there is no sign change within each band, as expected in 
our intuitive picture. 

In Fig[3Jb) , we fix —t' at 0.95^ and show the band 
filling dependence of T c , which takes its maximum around 
n = 1.25. This result shows that T c becomes low when 
n is too large, that is, when the Fermi level lies too far 
above the narrow band since this will make the system 
close to a purely single band model. T c also goes down 
when the Fermi level comes too close to or within the 
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narrow band B. 

In fact, when the Fermi level is within the narrow band, 
we find a t' dependence that is completely the opposite 
to what is seen in Fig^a). In Fig|3{c), we plot the 
eigenvalue A of the Eliashberg equation for n = 0.9 at 
T = 0.05i;. In this finite T c is not obtained for 

all of the t' values, so we plot A for a fixed temperature 
instead of T c . As seen in the figure, A takes its minimum 
around —if — ti- In this case, the Fermi level is within 
band B as can be seen from the Green's function shown 
in Fig|5Ja). Namely, the length between the two peaks 
is about 0.77T (which is in fact larger than the case of 
11 = 1.25 meaning that the band is not rigid), so that 
both of the bands has to be partially occupied^ We 
can see in FigOJa) that the Green's functions are small 
(the quasiparticle renormalization is strong) compared 
to those in Fig^a), so this should suppress T c . More- 
over, the gap function (Fig|5fc)) changes sign within the 
wide and the narrow bands near the Fermi surface, and 
this also should work destructive against superconductiv- 
ity. The sign change in the gap means that contributions 
from the pair scattering interactions within the narrow 
(and the wide) band are large, which is a consequence 
of the Fermi level crossing of the narrow band with a 
large density of states. This means that the peak struc- 
ture in the susceptibility (Figl^b)) originates from both 
interband and intraband scattering processes. 

The above results suggest that superconductivity with 
high T c may be obtained in the ladder compounds for the 
case of electron doping^ The condition for the maximum 
T c , —t' ~ i;, is of course unrealistic for the cuprate lad- 
ders, but we notice in Fig|3|that the enhancement of T c 
remains even if we deviate from —if = and relatively 
high T c (~ 0.04£;) is obtained even below — t' = 0.5^. 
Note that t' should be around t' ~ — 0.4i; in the ac- 
tual cuprate ladder compounds, assuming that t' has 
values similar to those for the layered cuprates such as 
YBaaCuaO^. 

To look in more detail into the possibility of high T c 
superconductivity for realistic values of the hopping in- 
tegrals, we plot in FigHJlthe band filling dependence of 
T c for —tf = 0.4t;. T c is not calculated near n = 1 be- 
cause (i) antiferromagnetic fluctuations strongly develop 
at high temperatures so that T c is not obtained, and in 
any case, (ii) FLEX loses its validity in the vicinity of 
half filling, where a Mott transition should take place. 
We find that a maximum T c of ~ 0.05t;, which is still 
considerably high, is reached at around n — 1.3, namely, 
when a large amount of electrons is doped. By contrast, 
for hole doping (n < 1), T c turns out to be much smaller, 
namely, of the order O.OOliz, which is of the same order 
as the T C (~12K) for the actual 14-24-41 compound^ The 
origin of the difference between the electron doped and 
the hole doped cases, can again be traced back to the 
Green's functions, the spin susceptibility and the gap 



functions, in which the characteristic features seen for 
if = — 0.95£/ (Fig^J still remain to some extent. The 
details on this point will be published elsewhere^ 

We have seen that our high T c mechanism due to co- 
existing narrow and wide bands works for a ladder sys- 
tem. Now, an interesting question is to ask how general 
this mechanism is. We have in fact found similar re- 
sults for the Hubbard model on a lattice where a pair of 
square lattices with in-plain nearest neighbor hoppings 
(t) is coupled by out-of-plain vertical (t v ) and diagonal 
(t') hoppings^ The Hamiltonian of this model is given 
as These results seem to suggest that the present mech- 
anism is likely to work on two band lattices having the 
Hamiltonian of the form 

fc,er 

where a and (3 are constants. Here, a flat band coexists 
with a wide band for a = ±1 (s(k) — — 2i;cos/c for the 
ladder, and — 2t(cosk x + cos fey) for the coupled planes). 
It would be an interesting future problem to investigate 
the validity of the present superconducting mechanism 
in a more wide class of models where wide and narrow 
bands coexist. 

To summarize, we have proposed a mechanism for high 
T c superconductivity in a two band system where wide 
and narrow bands coexist. When the Fermi level lies close 
to but not within the narrow band, T c becomes higher 
as the flatness of the narrow band increases, while com- 
pletely the opposite takes place when the Fermi level is 
within the narrow band. From this viewpoint, we have 
given a possible explanation for the rather low T c for the 
actual hole doped 14-24-41 ladder compound, and have 
further predicted a higher T c for the case of electron dop- 
ing. As for the 14-24-41 compound, recent NMR experi- 
ments have observed a coherence peak and an unchanged 
Knight shift across T c ~ As for the presence of the co- 
herence peak, the singlet gap function obtained in our 
study can be consistent with the experiments since the 
gap does not change its sign on the Fermi surface. As for 
the unchanged Knight shift, we believe further investiga- 
tion is necessary, but even this is due to spin-triplet pair- 
ing originating from effects not included in the present 
study (e.g., phonons), we have to explain why T c of the 
singlet pairing is even lower than 12K. The present view- 
point can still be relevant for answering this question, 
and we believe that the prediction for a higher T c in the 
case of electron doping remains valid. 
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